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Abstract 

We extend the study of velocity quantization phenomena recently found in the classi- 
cal motion of an idealized ID model solid lubricant - consisting of a harmonic chain 
interposed between two periodic sliding potentials [Phys. Rev. Lett. 97, 056101 
(2006)]. This quantization is due to one slider rigidly dragging the commensurate 
lattice of kinks that the chain forms with the other slider. In this follow-up work 
we consider finite-size chains rather than infinite chains. The finite size (i) permits 
the development of robust velocity plateaus as a function of the lubricant stiffness, 
and (ii) allows an overall chain-length re-adjustment which spontaneously promotes 
single-particle periodic oscillations. These periodic oscillations replace the quasi- 
periodic motion produced by general incommensurate periods of the sliders and the 
lubricant in the infinite-size model. Possible consequences of these results for some 
real systems are discussed. 



1 Introduction 



The present paper extends the study of a one- dimensional (ID) non-linear 
model, inspired by the tribological problem of two sliding surfaces with a thin 
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solid lubricant layer in between [1], to the case of a lubricant ID "island" of 
finite size. Previous work [1,2] found robust, universal and exactly quantized 
asymmetric velocity plateaus in the classical dynamics of an infinite-size chain 
subject to two relatively sliding periodic potentials. The infinite chain size was 
managed - in the general incommensurate case - by means of periodic bound- 
ary conditions (PBC) and finite-size scaling. The plateaus of chain velocity 
as a function of several model parameters were shown to be due to motion 
of kinks (topological solitons), which generally exist in a chain submitted to 
a periodic potential. These nonlinear excitations generated by the first, sta- 
tionary periodic potential are set into motion by the external driving which is 
provided by the second periodic potential, sliding with velocity v^xt- While the 
chain kinks are thus dragged with velocity fcxtj the overall chain velocity is 
smaller, and fixed by the kinks nature and density. That in turn depends only 
on the ratio of the period of one slider to that of the chain, which is dictated 
by the interparticle spacing, and enforced by the PBC. 

The present work considers a finite open-boundary chain, such as for example 
a hydrocarbon chain, or a graphite fiake interposed between two sliding crystal 
faces [3]. Unlike the infinite chain, the open-boundary chain can elongate or 
shorten, at the cost of some harmonic potential energy, effectively modifying 
its linear density, and thus the kink density which is the relevant length ratio. 
We find that it does indeed elongate or shorten so as to realize a precise 
commensurate relation to the other slider. This adaptive relaxation is such 
as to produce perfectly periodic oscillations of the single particles superposed 
with their average drift at the quantized velocity. 



2 The model 



Using the same language and notation of previously studied confined models 
[4,5,6,7,8], we represent a solid lubricant layer as a chain of harmonically 
interacting particles interposed between two rigid generally (but not necessar- 
ily) incommensurate sinusoidal substrates (the two '"sliding crystals" , sketched 
in Fig. 1) externally driven at a constant relative velocity Vext- The equation 
of motion of the i-th particle is: 



mxi = — - [F+ sin k^{xi — v^t) + F_ sin /c_(xi — v_t)] 

+ K{xi+i + Xi_i-2xi)--iY^{xi-v±) , (1) 

± 

where m is the mass of the particles, K is the chain spring constant, and 
k± = 271 /a± are the wave- vector periodicities of potentials representing the 
two sliders, moving at velocities v±. We set, in full generality, u+ = and 
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V- — Vext- 7 is a phenomcnological parameter substituting for various sources 
of dissipation, required to achieve a stationary state, but otherwise playing 
no major role in the following. F± are the force amplitudes representing the 
sinusoidal corrugation of the two shders (we will commonly use F-/-F+ = 1 
but we checked that our results are more general). We take a+ — 1, m — 1, 
and = 1 as our basic units, and all quantities are measured in suitable 
combinations thereof. The relevant length ratios [9,10] are defined by r± = 
a±/ao; we assume, without loss of generality, r_ > r+. The inter-particle 
equilibrium length ag enters explicitly the equations of motion (1) of the first 
(i = 1) and last {i — N) particle whose restoring force terms in Eq. (1) are 
K{x2 — xi — ao) and K{xn-i — xn + oo), respectively; this implements open 
boundary conditions (OBC). 

Upon sliding the substrates, Vext 7^ 0, the lubricant chain shdes too. Despite 
the apparent generic symmetry between the two sliders, the time-averaged 
chain velocity w = t'cm/'^cxt) is generally asymmetric, namely different from 
1/2. In a previous study on this model [1,2] it was shown that, for an infinite 
chain and periodic boundary conditions (PBC), w is exactly quantized, for 
large parameter intervals, to plateau values that depend solely on the chosen 
commensurability ratios (r+,r_). As the present finite-size OBC simulations 
will show, the PBCs are not crucial to the plateau quantization, which occurs 
even for a lubricant of finite and not particularly large size N. The main 
difference between OBC and PBC is that, while in PBC the chain length, 
and thus the length ratios r±, are fixed, in OBC the chain can lengthen or 
shorten with respect to its equilibrium size. We find that during sliding the 
chain length gradually reaches a natural attractor value and oscillates around 
it. We define a new effective inter-particle length and corresponding length 
ratio 



where L = {xn — xi) is the average chain length after the initial transient. 
The effective length ratio plays a central role in the understanding of the 
velocity plateaus of the finite-size lubrication model. 



3 Results and theory 



The driven dynamics of the lubricant is studied by integrating the equations of 

motion (1) starting from fully relaxed springs {xi = iaQ, Xi = fcxt/2), using a 
standard fourth-order Runge-Kutta method. After an initial transient where 
length relaxation takes place, the system reaches its dynamical stationary 
state, at least so long as 7 is not exactly zero. Figure 2 shows the resulting 
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time-averaged center-of-mass (CM) velocity function of the chain 

stiffness K, for an irrational choice of (r+, r_), and two values of A^, defining a 
relatively short chain [N = 15) and one of intermediate length (A^ = 100). We 
find that w is generally a complicated function of K, with flat plateaus and 
regimes of continuous evolution, not unlike what is found in the inflnite-size 
limit studied through PBC simulations [1] . 

To investigate the peculiarities brought about by finite size, we analyze the 
dynamics for a large number of values of (r+,r_) and K, and observe that: 
(i) one or more velocity plateaus as a function of K occur for large ranges of 
(r+,r_); (ii) the velocity ratio w of most plateaus satisfies 



This result should be compared with the relation (w = 1 — r^^) vahd for 
the main plateau of PBC calculations [1], where the length ratios (r+,r_) 
are fixed: in Eq. (3) the new effective length ratio r"^ of Eq. (2) replaces r+. 
Figure 3 collects the observed plateau velocity ratios for a range of values of 
r_|_ and for fixed r_ = ^/lOl. Clear trends emerge: plateau velocities depend 
continuously on r+, thus indicating that also the effective chain length ratio 
evolves continuously with r+ in finite ranges. Several plateaus appear to 
follow different curves. 

These data arc conveniently organized and understood by plotting r'^ rather 
than w, as a function of r+, as is done in Fig. 4. ^ All points fit perfect straight 
lines through (r_|_,r^^) = (0,1). Most lines have slope g/r_ with integer q. 
Occasional plateaus fit this relation with half-integer q (—5/2, 7/2, and 9/2, 
in the calculations of Fig. 4) . 

Several calculations carried out with different values of r_ confirm that in 
general the ratio satisfies 

rf = l + g!± (4) 



with q taking simple fraction (often integer) values. This general behavior 
indicates that the plateau dynamics leads the finite-size lubricant toward a 
dynamical configuration where not only its velocity but also its length is quan- 
tized. 



For reasons of numerical stability, this figure reports r'^ obtained by inversion 
of Eq. (3), but the same plot could have been obtained directly by applying the 
definition (2). This equivalence is confirmed by the comparison of the two methods 
shown in Fig. 4 for 2 points. 
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We can understand this phenomenology in terms of kinks (i.e., local compres- 
sions of the chain with substrate potential minima holding more than just one 
particle [11]), as described in [1]. Assume initially integer q. The basic hypoth- 
esis explaining the relation (4) are: (i) the particles tend to singly occupy the 
a+-spaced minima of the bottom potential, with occasional kinks to release the 
spring tension; (ii) kinks group in bunches, each sitting in a period a_ of the 
top substrate; (iii) kink bunches are g-fold, i.e. they collect q individual kinks 
(negative q indicates the number of anti- kinks). After the initial transient the 
chain length becomes on average very close to L — {N — N^ink — 1) a+. The 
number of kinks thus equals the number q of kinks per bunch times the total 
number L/a_ of bunches in the chain: Nkink — qL/a-. By eliminating Nkink, 
we obtain 

a- + qa+ 



This is consistent with an average inter-particle distance 

eff L ^- 1 lci\ 



and thus with the effective length ratio of Eq. (4) . In general, for non- integer 
Q = 'nk/n^ values, this interpretation remains valid: bunches of a total of 
rik kinks distribute themselves in n_ minima of the a_ lattice, q indicates 
therefore the density (coverage fraction) of kinks on the a_ lattice. 

The dynamically stable plateau attractors of the open-boundary chain are 
therefore characterized by a lattice of kinks perfectly commensurate to the a_ 
lattice. In the infinite-size PBC model, a rational kink coverage corresponds 
to commensurate encounter frequencies /+ and /_ of the generic lubricant 
particle with the two substrates, which in turn occurs for very special values 
of r± [1]: these (r+,r_) values are characterized by perfectly periodic single- 
particle dynamics. It is rather remarkable that the open-chain model, without 
the necessity of any careful fine tuning of (r_|_ , r_ ) realizes a self-organized 
commensurate kink lattice automatically producing perfectly periodic single- 
particle oscillations in a generally incommensurate context. 

Note that the points of Fig. 4 lying along the g = line indicate perfect 
commensuration between the lubricant chain and the bottom substrate, to 
which the chain remains pinned: this can be realized by paying a moderate 
harmonic-energy cost only close to r+ = 1. For the same reason, all plateaus 
appear close to the r+ = line (dot-dashed hne in Fig. 4). 

For the parameters of Fig. 3, nontrivial plateaus are found only for 0.7 < 
r+ < 1.6. Outside this range, relation (4) produces r^^ , for small integer q, 
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very much different from r+. Smaller values of r_ ~ 2 (rather r_ ~ 10 as in 
Fig. 3) generates an analogous set of plateaus for r+ > 1.6. 



4 Size dependence 

As Fig. 2 suggests, strong size effects are observed, especially for large K. In 
particular, the "natural" symmetric large- limit w = 1/2 found with PBC 
in Ref. [1] is rarely reached using OBC: for N — 100 the chain is pinned to the 
bottom substrate {vcra — 0), while for = 15 it moves at Vcra — 't'ext- Figure 5 
shows the dependence of the velocity ratio w on the particle number A^, for 
a fixed length ratio r_|_ and two different values of r_, in the large- A' limit. 
We observe that the chain is pinned to the a+ substrate when A^ occurs to be 
(nearly) a multiple of the length ratio r_. In all other cases, the chain follows 
the a_ substrate, at velocity T;ext- 

This changing behavior can be understood as follows. For large spring stiffness 
the kink dynamics is suppressed, the lubricant particles placing themselves 
at nearly regular distances ~ Cq. It is energetically favorable for a chain shorter 
than a_ to sit in a minimum of the top potential and then stick to it. Even 
if the chain is longer than a_, its length is generally not an exact multiple of 
the periodicity of the top substrate. For this reason, a finite end part of the 
chain remains similarly trapped in the minima of the top substrate: this is 
what occurs for most A^ values in Fig. 5. On the other hand, when A"ao is a 
multiple of the a_ period (i.e. when A" is close to a multiple of r_ = a_/ao), 
minima and maxima of the upper potential compensate each other, so that 
there is no preferential relative position of chain and top substrate. For such 
special sizes, the chain remains weakly pinned to the bottom substrate, as 
illustrated in Fig. 5. The values A" = 15 and A" = 100 of Fig. 2 represent the 
two situations. 

Other intermediate values occur for specific sizes, but the finite-size scaling 
for large A" is obviously non-trivial. While Fig. 2 indicates that for small and 
moderate K, size effects, if any, are very small, at large K they affect the 
dynamics substantially. Further work is needed to understand the large-X 
size-scaling in detail. 



5 Discussion and Conclusions 

We have shown that chains of finite and even small size, driven in between two 
periodic sliders, move with characteristic quantized velocities, much like the 
infinite-size ones do. We find that a finite-size chain length re-adjusts in such 
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a way as to realize a lattice of topological solitons (kinks) commensurate to 
the period of the smoother slider. A consequence of this sclf-commensuration 
is a periodic single-particle motion, even for incommensurate initial choices of 
the periods. 

The likely reason behind this phenomenology is as follows. Consider initially a 
single periodic slider of length ratio r+ and periodic boundary conditions. The 
lubricant chain will then form a regular lattice of kinks, which repel each other. 
When the second slider is introduced, the lattice of kinks will be generally 
incommensurate with the period a_, and this brings an irregular distribution 
of bunches of kinks, which on average reconstruct the correct density of kinks. 
An OBC finite chain is able to relax, by paying some extra harmonic strain, 
in such a way as to enforce an optimal local density of kinks which satisfies 
both periodic sliders. 

The phenomena just described for a model ID system are unique, and it would 
be interesting if they could be observed in real systems. Nested carbon nan- 
otubes [12], or confined one-dimensional nanomechanical systems [13], are one 
possible arena for the phenomena described. Though speculative at this stage, 
one obvious question is what aspects of the phenomenology just described 
might survive in two- dimensions (2D), where tribological realizations, such as 
the sliding of two hard crystalline faces with, e.g., an interposed graphite flake, 
are conceivable. Our results suggests that the lattice of discommensurations 
- a Moire pattern- formed by the fiake on a substrate, could be dragged by 
the other sliding crystal face, in such a manner that the speed of the flake 
as a whole would be smaller, and quantized. This would amount to the slider 
"ironing" the kinks onward. Dienwiebel et al. [3] demonstrated how incom- 
mensurability may lead to virtually friction-free sliding in such a case, but no 
measure was obtained for the flake relative sliding velocity. Real substrates 
are, unlike our model, not rigid, subject to thermal expansion, etc. Never- 
theless the ubiquity of plateaus shown in Fig. 2, and their topological origin, 
suggests that these effects would not remove the phenomenon. 
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Fig. 1. (Color online) A sketch of the model with the two periodic sliders (periods 
a+ and a_) and the lubricant harmonic chain of rest nearest-neighbor length oq- 
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Fig. 2. (Color online) Average drift velocity ratio w = Ucm/fext 
as a function of the lubricant chain spring stiffness K for 
(r+,r_) = (^51+ J^,VWI) ~ (1.2249,10.050) and for two different finite 
sizes of the chain = 15 and 100. OBC are applied, and 7 = 0.1, Vext = 0.01. The 
plateau labeling refers to values of q, defined in Eq. (4). Note the similar behaviour 
for if < 10 in spite of the different chain size, and the different limit for large K, 
due to finite-size effects. 
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Fig. 3. (Color online) The plateau velocity ratio if as a function of r_|_, for 
r__ = 101^''^. Circles, diamonds and triangles represent respectively the first, second 
and third plateau found for increasing \w\. For comparison, we report the plateau 
velocity ratio w = 1 — r^^ (dot-dashed curve) of the fixed- length PBC calculations 
[1]. Here 7 = 0.1, Uext = 0.01, and a finite chain of = 100 particles (OBC) is 
considered. 
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Fig. 4. (Color online) The effective length ratio r'^^ = a+/aQ^, calculated from the 
velocity data of Fig. 3 using = (1 — w)~^ , the inverse of Eq. (3). The stars are 

values obtained based on the definition (4) by time-averaging the chain length, 
showing perfect agreement with the velocity-derived values. The dot-dashed line 
represents the identity = (undeformed chain length). Solid straight lines: 

= 1 -|- g for integer q. 
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Fig. 5. (Color online) Size scaling of the average drift velocity ratio w for 
= ^.^^ = (f), and two different values of the ratio r_ = (101)^/^ ~ 10.05 and 
r„ = (101)^/^ + 3 ~ 13.05. The number of particles N for which the chain is pinned 
to the a+ substrate {w = 0) are close to integer multiples of r_. All calculations 
refer to 7 = 0.1, Vext = 0.01, and K = 1000, a representative of the large- -ftT limit. 
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